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ON  ELECTRON  BEAM  IONIZATION  OF  AIR  AND  CHEMICAL  REACTIONS  FOR 

DISTURBED  AIR  DEIONIZATION 

1.  INTRODUCTION 

When  a  pulsed  high  energy  electron  beam  propagates  in  the  atmosphere, 

It  loses  energy  through  several  physical  processes.  The  energy  lost  by  the 
beam  is  absorbed  by  the  air  molecules  resulting  in  a  certain  degree  of 
heating  and  ionization  of  the  atmosphere.  The  degree  of  heating  and 
ionzation  depends  on  the  amount  of  energy  deposited  in  the  atmosphere  by  the 
electron  beam. 

The  heated  and  the  partially  ionized  atmosphere  (disturbed  air)  will 
tend  to  relax  to  its  ambient  condition  after  the  passage  of  the  electron 
beam.  The  relaxation  and  the  deionization  of  a  disturbed  air  proceed  through 
a  host  of  atomic,  molecular  and  chemical  processes. 

This  report  will  deal  with  the  generation  of  the  disturbed  air  and  its 
subsequent  relaxation.  It  describes  methods  for  detailed  treatment  of  the 
energy  deposition,  the  deionization  processes  and  the  air  chemistry  to 
describe  the  conductivity  of  an  ionized  air.  i'he  report  also  treats  the  role 
of  water  vapor  on  the  air  chemistry  and  reviews  in  detail  the  current  status 
of  the  reaction  rates  pertinent  to  the  modeling  and  the  understanding  of  the 
deionization  phenomena  in  air. 

The  electron  energy  deposition  in  air  clearly  will  generate  a  plasma 
where  a  large  number  of  plasma  processes  will  also  occur.  However,  this 
report  will  deal  with  the  electron  beam  interaction  with  air  purely  from  the 
point  of  view  of  the  basic  atomic,  molecular  and  chemical  reactions  and  their 
role  on  the  ionization  and  deionization  in  general. 

2.  ELECTRON  BEAM  GENERATED  IONIZATION  IN  AIR 

We  consider  high  energy  (E  ^  1  MeV)  and  high  current  (I  *  lkA)  pulsed 
electron  beams  propagating  in  the  atmosphere.  These  beam  electrons  lose 
energy  by  three  distinct  processes  which  result  in  the  ionization  of  the 
air  species.  The  energy  loss  processes  are;  the  inelastic  collisions,  the 
Bremsstralung  radiation  and  the  ohmic  heating,  which  require  distinct  energy 
deposition  treatments  to  calculate  their  contributions  to  the  total  ionization 
in  air. 

In  this  section,  the  electron  beam  energy  loss  processes  will  be 
discussed  briefly  along  with  the  energy  deposition  methods  which  can  be  used 
to  calculate  the  ionization. 
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2.1  ENERGY  LOSS  BY  COLLISIONS 

The  primary  electrons  (beam  electrons)  lose  energy  by  Inelastic 
collisions  with  air  molecules.  The  lost  energy  is  expended  in  excitations , 
dissociations,  ionizations  of  the  air  specie')  and  the  generation  of  the 
secondary  electrons.  The  energetic  secondary  electrons  undergo  the  same 
energy  loss  processes  and  hence  contribute  to  the  total  collisional 
ionization. 

The  energy  loss  per  unit  path,  dE/dx,  by  relativistic  electrons  due  to 
collisions  with  neutral  species  of  density  N  can  be  expressed  as1 


where  E  is  the  total  energy  of  the  electron,  T  is  its'  kinetic  energy,  @  -  v/c, 
v  is  the  velocity  of  the  electron,  Z  is  the  nuclear  charge  of  the  scattering 
atom  and  I  is  the  mean  excitation  energy  of  the  atom,  while  the  res :  of  the 
symbols  have  their  usual  meanings.  Equation  (1)  can  be  expressed  as: 


4.93  x  10“25  me2  N  Z 

-  (2Y  -  1)  log  2  + 


(y  l)2  (y  +  1)  (me2)2 
2I2 

(2) 


where  y  is  the  relativistic  factor  (E  ■  ytuc2).  For  relativistic  electrons 
in  air,  at  the  sea  level,  the  energy  loss  per  unit  path  in  units  of 


eV,  is 


2 


dE 

dx 


97.33 

|c  ,2 


V*  -  1 


log  | 1.47  x  107  (v  -  l)2  (y  +  i)j 


(2V  -  1) 


1  (Y  -  1)! 


log  2  +  —  + 

.  ,2 


(3) 


where  Z  «  7*2  and  1-94  eV1,s  were  utilized  in  Equation  (2).  For  highly 
relativistic  electrons  the  first  term  in  the  bracket  in  Equations  (1)  to  (3) 
is  the  dominant  term  and  it  suffices  to  use  this  term  alone  for  the  purposes 
of  calculations.  In  Figure  1  we  show  the  energy  loss  due  to  collisions  by 
electrons  in  air  at  the  sea  level  using  Equation  (3).  However,  it  should  be 
noted  that  Equations  (1)  -  (3)  do  not  Include  the  density  effect  or  what  is 
generally  called  the  polarization  effect.  The  inclusion  of  this  effect  which 
could  be  calculated,3  would  result  in  reducing  the  energy  loss  per  unit  path* 
Obviously,  the  polarization  effect  becomes  appreciable  starting  with  electrons 
with  energy  of  20  MeV  and  higher.  The  energy  loss  by  collision  with  polariza¬ 
tion  effect  included,  is  also  shown  in  Figure  1  where  data  for  E  s  100  MeV 
.and  lower  are  taken  from  Reference  4. 

Some  of  the  features  of  the  energy  loss  through  collisions  can  be  seen 
from  Figure  1,  where  the  energy  loss  decreases  sharply  with  increasing  energy 
until  it  reaches  a  minimum  where  E  2  MeV,  then  starts  to  rise,  but  slowly 
(logarithmically),  with  increasing  energy.  The  polarization  correction  makes 
the  etvergy  loss  by  collisions  almost  a  constant  beyond  E  ^  200  MeV. 

2.2  ENERGY  '..OSS  BY  RADIATION 

High  energ>  electrons  loss  their  energy  also  by  radiation  as  they 
traverse  a  gaseous  medium  or  a  material  target.  The  radiative  energy  loss  is 
due  to  the  interaction  of  the  Incident  electron  with  the  Coloumb  fields  of  the 
nucleus  and  the  atomic  electrons.  The  acceleration  and  deceleration  effects 
on  the  incident  ele:tron  result  in  Bramsstrahlung  radiation.  The  electron 
energy  loss  by  radiation,  in  air,  becomes  appreciable  when  E  5  40  MeV  and 
rises  sharply  with  increasing  electron  energy  (see  Figure  1),  and 
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KINETIC  ENERGY  (MeV) 


becomes  the  dominant  energy  loss  process  for  highly  relativistic  electrons 
when  a  200  MeV. 

The  energy  loss  by  radiation  for  a  relativistic  electron  with  energy 

E  is1 

8  r  -  N  J  hv  <t>  (E,  s>)  dv  (4) 

where  v  is  the  frequency  of  the  emitted  radiation  and  ve  is  the  upper  limit 
for  the  radiated  frequency  defined  as  E  -  hv^  +  me  2  and  <t>  (E,  v)  dv  is  the 
cross  section  for  the  production  of  a  photon  with  frequency  v  in  the  range  of 
v  +  dv.  These  cross  sections  have  been  calculated  by  Bethe  and  Heitler  (see 
reference  1)  and  are  given  over  a  wide  range  of  the  screening  parameter.1*6 
As  an  example  for  complete  screening1*6  (extremely  relativistic  electrons), 
the  Bremsstrahlung  radiation  production  cross  section  is1*6 


where  r0  is  the  electron  radius.  Using  Equation  (5)  into  Equation  (4)  and 
performing  the  Integration,  one  obtains 
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which  can  be  written  as 


where,  the  radiation  length,  L  ,  is  defined  as 

r 


(7) 


l“l  ■  5*7  x  10-“  (z2  +  Z)  N 


(8) 


For  highly  relativistic  electrons,  the  radiation  length  in  air,  at  the  sea 

level,  is  300  meters.  This  implies  that  the  highly  relativistic  electrons 

traversing  300  meters  of  air  at  the  sea  level  will  have  their  energies 

reduced  to  1/e  of  their  original  value;  the  loss  being  due  entirely  to 

radiative  emission,  i.e.  Bremastrahlung.  The  energy  loss  per  unit  length 

(eV/cm)  by  the  relativistic  electrons  in  air,  at  the  sea  level  is  shown  in 

Figure  1,  where  Equation  (6)  is  utilized  with  E  *  100  MeV.  For  E  <  100  MeV, 
dE  \ 

the  ^  j  values  are  taken  from  Reference 
2.3  ENERGY  LOSS  BY  OHMIC  HEATING 

In  addition  to  the  energy  loss  by  inelastic  collisions  and  Bremsstrahlung 
radiation,  the  relativistic  beam  electrons  lose  energy  through  the  ohmic 
heating  of  the  plasma  electrons.  The  plasma  electrons  gain  energy  from  the 
beam  induced  electric  field  which  drives  a  return  current.  The  enhancement 
in  the  plasma  electron  energy  leads  to  the  avalanche  ionization7  of  air, 
especially  under  reduced  air  density  conditions.  When  the  air  density  is 
low  (few  Torrs  and  lower)  additional  processes  such  as  the  two  stream 
instability, 8,9  will  contribute  to  the  generation  of  the  plasma  electrons 
and  hence  act  a &  another  mechanism  for  the  energy  loss  by  the  beam. 

The  generation  of  the  plasma  electrons  by  collisional  and  other  means 
results  in  an  air  plasma  with  conductivity  ?  given  by 


where  v  Is  the  momentum  transfer  collision  frequency  which  includes 

electron-neutral  and  electron- ion  collisions. 

Whan  the  beam  radial  electric  field  is  much  less10*11  than  the  axial 

electric  field,  E  .  the  return  current  density,  J  ,  is  related  to  E  by 
z  z  z 


j 


z 


a  E 

z 


(10) 


To  drive  the  return  current,  the  beam  electrons  lose  energy  by  ohmic  heating. 
The  energy  lost  per  unit  length,  assuming  a  uniform  current  density,  is 


dE 

dx 


ohm 


t 


A 


/j  E  dt 
z  z 


(ID 


where  A  is  the  beam  area  cross  section.  The  parameter?  a,  j  and  E  in 

z  z 

Equation  (11)  can  be  obtained  using  the  appropriate  rate  equations  to 
calculate  ct  and  the  usual  circuit  equations10  for  the  current  driveu  by  the 
electron  beams.  However,  an  estimate  can  be  made  for  the  energy  loss,  due 
to  ohmic  heating.  Using  a  one  dimension*!  circuit  equation,  assuming  a  fast 
current  rise  time  and  beam  pulses  short  enough  so  that  Ir  is  roughly  constant, 
one  obtains, 


dE 

dx 


ohm 


I 

n 


(12) 


where  L  is  the  circuit  inductance  per  unit  Ipngth,  1^  and  1^  are  the  beam 
and  net  currents,  respectively. 

2.4  IONIZATION  CALCULATION 

In  this  section  we  discuss  methods  for  the  ionization  calculations  due 
to  the  energy  loss  processes  presented  in  Sections  2.1  -  2.3. 

The  collisions!  ionization  consists  of  the  direct  ionization  by  the 
beam  electrons  (primary  electrons)  and  the  ionization  by  the  secondary 
electrons  generated  by  the  act  of  the  primary  ionization.  One  of  the  often 

7 


idui&uiai 


used  concepts  £or  the  collisional  Ionization  calculation  is  the  method  of  the 
continuous  slowing  down  approximation.12"14  In  this  method  one  calculates 
the  primary  ionization,  obtains  the  secondary  and  the  tertiary  electron 
distribution  and  their  contributions  to  the  total  collsional  ionization.  It 
has  been  utilized14”13  for  electron  and  proton  depositions  in  Na»  02  and  0 
in  the  upper  atmosphere  for  the  understanding  of  the  auroral  phenomenon.  This 
energy  deposition  concept  has  been  improved  to  include  discrete  energy 
depositions  for  electrons  with  energy  below  500  eV  and  has  been  applied  to 
auroral  phenomenon17*18  as  well  as  depositions  in  other  gaseous  elements.13"21 
Such  a  treatment  of  a  detailed  energy  deposition  is  in  progress  At  NRL.  It 
requires  a  large  number  of  ionization  and  excitation  cross  sections  and  has 
the  advantage  of  providing  excitation  rates  to  all  electronic  states,  for  a 
detailed  air  chemistry  and  plasma  emission  calculations.  Some  of  these 
electronic  states,  with  radiation  emissions  in  the  optical  region,  can  be 
utilized  for  various  beam  diagnostic  purposes. 

The  Ionization  of  air  by  the  Bremastrahlung  radiation  is  more  complicated 
and  depends  on  the  flux  and  the  energy  distribution  of  the  emitted  photons. 
These  photons,  absorbed  and  scattered  by  the  air  molecules,  generate  free 
electrons  by  (a)  ,-hotoabsorption,  (b)  Compton  effect  and  (c)  pair  production. 
The  free  electrons  generated  in  this  manner  deposit  their  energies 
collisionally  as  well  as  by  radiative  losses.  Thus,  the  Bremsstrahlung 
iouizatlon  of  air  can  be  described  and  calculated  as  a  cascade  shower5*22 
ionization. 

A  detailed  Bremastrahlung  energy  deposition  in  air  is  accomplished  by 
the  calculation  of  the  flux  and  the  energy  distribution  of  the  emitted 
photons  for  each  primary  electron.  These  photons  then  are  deposited  in  air 
as  a  function  of  distance  using  the  cross  sections  for  the  ionization  of  air 
due  to  the  photoabsorption,  Compton  effect  and  pair  production.  These  cross 
sections,  shown  in  Figure  2  in  the  form  of  the  attenuation  coefficient,23  11, 

(y,  -  y.^,  y,^  »  N^)  indicate  the  role  or  each  process,  in  air,  as  a 

function  of  the  radiation  wavelength.  For  each  process,  the  energy  of  the 
ejected  electron  is  calculated  and  its  energy  is  deposited  according  to  the 
energy  loss  processes.  Such  an  energy  deposition  for  the  Bremsstrchlung 
radiation  is  under  consideration  at  NRL. 

The  calculation  of  the  beam  energy  loss  by  ohmic  heating  can  be  carried 
out  also  ir.  detail  by  solving  Maxwell's  eq'iations,  the  kinetic  equations  and 
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the  appropriate  circuit  equation  to  calculate  various  currents  (Iq,  1^,  etc.) 
and  the  air  conductivity. 

These  energy  deposition  schemes  described,  briefly,  to  account  for  the 
ionisation  in  air  were  described  separately.  In  reality,  however,  the  energy 
loss  processes  are  simultaneous  occurrences,  albeit  oue  more  important  than  the 
other  at  different  tines,  pressures,  and  positions  of  interest.  Furthermore, 
good  estimates  for  ionization  can  be  obtained  using  simpler  methods  e.g.  the 
energy  per  ion  pair  generation  concept  for  the  collisional  and  radiative 
losses  and  using  a  simple  circuit  equation1'1  to  describe  the  ohmic  heating. 

3.0  THE  I0MIZATI0N  AND  THE  DEIONIZATION  PROCESSES  IN  AIR 

In  this  section,  we  give  an  overview  of  the  ionization  and  the 
deionization  processes,  in  air,  which  occur  during  the  passage  of  a  pulsed 
relativistic  electron  beam.  The  purpose  of  this  overview  is  to  provide  the 
understanding  of  the  patterns  of  the  charge  flow  and  how  they  control  the 
time  history  of  the  electron  density  and  consequently  the  air  conductivity. 

Tne  ionization  in  Che  ambient  air  proceeds  by  electron  impact  on  the  air 
molecules  and  by  their  absorption  of  radiation.  The  collisional  ionization 

4*  ^  ^  *4* 

produces  air  species  ions  in  the  form  of  Na  ,  H  ,  02  and  0  .  These  ions  are 
formed  by  two  distinct  processes,  which  are  the  ionization  and  the 
dissociative  ionization  of  the  aic  molecules,  as  described  by  Equations  13- lo. 


e  +  N2 

N2+  +  2e 

(13) 

e  +  Na 

N+  +  N  + 

2e 

(lb) 

e  +  02 

■4 

02+  +  2e 

(15) 

e  +  02 

-♦ 

0+  +  0  + 

2e 

(16) 

The  ionization  resulting  from  the  absorption  of  radiation  proceeds  by 
photoionization  and  by  collisional  Ionization,  where  the  last  process  is  due 
to  the  Compton  and  the  pair  produced  electrons.  For  photons  with  low  energies 
(hv  a*  15  -  200  eV) ,  the  ionization  proceeds  by  the  ionization  of  the  valence 
electrons  as  described  la  Equations  17  -  20 

hv  +  Na  ^  N2"*"  +  e  (17) 
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(iQ) 


hv  +  N2  -»  N+  +  M  +  e 
hv  +  Oa  -♦  02^  +  e  (19) 

hv  +  Oa  -4  0+  +  0  +  e  (20) 

However,  for  photons  with  energies  above  200  eV  and  where  Compton  effect  is 
small  the  photoicnization  proceeds  via  the  inner  shell  ionization  and  one 

|  |  —a 

generates  N2  and  03  due  to  the  Auger  process24  because  the  x-ray 
fluorescence  yield25  is  very  low  for  low  z  elements.  These  doubly  ionized 
molecules,  however,  will  be  converted  to  2  N2  and  2  O2  through  charge 
exchange  with  their  respective  molecules. 

The  oxygen  and  nitrogen  atoms  are  also  generated  by  the  electron 
Impact  dissociations  of  N2  and  O2 

e  +  N2  -*  N  +  N  (21) 

e  +  Oa  «*  0+0  (22) 


and  by  the  dissociative  recombinations  of  the  molecular  ions 


e  +  N2+  N  +  N 

(23) 

e  +  O2  0+0 

(2k) 

e  +  N0f  -♦  N  +  0 

<25) 

Since  these  processes  (i.e.  Reactions  -1-25)  occur  during  the  electron  beam 
pulse,  the  atomic  species  are  ionized  by  electron  impact  to  also  produce 
the  ionic  species  of  N  and  0  .  The  *  .tie  oxide,  NO,  is  produced  mainly  by 

N  +  02  h  NO  +  0 

(26) 

N(®D)  +  Ls  ■*  NO  +  0 

(27) 

The  N0+  ion,  however,  is  produced  by  electron  Impact  and  photoionization  of 
NO,  by  charge  exchange  and  by  ion-molecule  rearrangement  processes  which  are 


11 


described  below 


(28a) 


e  +  NO  -♦  NO**"  +  2e 
hv  +  NO  -»  NO+  +  e 


^  1 

* 

N ,  0+,  0a+ 
l  J 

+  NO  -»  N0+  f 

N,  0,  02 

N+  +  Oa  NO+  +  0 
0+  +  Na  -»  NO"*"  +  N 


(28b) 

(29) 

(30) 

(31) 


The  deionization  of:  the  air  plasma  occurs  during  the  passage  of  the  electron 
beam  and  continues  afterwards,  until  the  free  electrons,  the  positive  ions 
and  the  negative  ions  have  disappeared  completely. 

The  deionization  proceeds  as  follows;  the  electrons  disappear  in  two  ways 

(1)  by  recombinations  with  the  ions,  (2)  by  attachment  and  dissociative 

attachments  to  molecules.  The  attachments  generate  negative  ions  which 

deplete  the  positive  ion  concentrations  by  mutual  neutralization  (ion-ion 

recombination).  However,  during  the  deionization  the  positive  ions  undergo  a 

multitude  of  reactions,  such  as  charge  exchange,  ion-molecule  rearrangements 

and  associations.  These  processes  are  Important  because  generally  they 

convert  an  atomic  ion  into  a  molecular  one  or  a  molecular  ion  into  a  cluster 

ion,  thereby  speeding  up  the  plasma  recombinations.  This  is  because  the 

molecular  and  cluster  ions  recombine  faster  with  electrons  compared  to  the 

recombination  of  electrons  with  atomic  ions,  especially  for  air  plasmas  with 

N  <  4  x  1016  cm”? 
e 

To  illustrate  the  deionization  processes  in  air,  the  molecular  ions 
N2  and  02  recombine  disoociatively  with  the  plasma  electrons  i.e. 

e  +  Na+  -»  N  +  N  (23) 

e  +  0a+  -»  0  +  N  (24) 


These  ions,  also  undergo  collisional-radiative  recombination  along  with  atomic 
ions.  However,  for  regions  of  interest  we  can  Ignore  the  radiative 
recombination  and  consider  only  the  three-body  (collisional)  recombination, 
e.g. 

e  +  e  +  (N+,  Na+)  -♦  e  +  (N,  Na)  (32) 


12 


(33) 


e  +  a  +  (0  ,  Og  )  e  +  (0,  Og) 

Three-body  recombination  also  proceeds  with  the  third  body  being  a  neutral 
species,  e.g«, 

e  +  M  +  (N+,  Ng+)  ->  M  +  (N,  Ng) 

e  +  M  +  (0+,  0g+)  >4  M  +  (0,  Og) 

Apart  from  these  direct  paths  which  deplete  the  free  electron 
concentration,  the  free  electrons  and  the  atmospheric  ions  follow  different 
paths.  The  free  electrons  are  lost  by  three-body  and  dissociative  attachments 
with  Og,  whereupon  stable  negative  ions  are  generated,  e.g., 

e  +  02  +  M  *  0a“  +  M  +  (36) 

e  +  Og  0“  +  0  (37) 


(34) 

(35) 


The  positive  ions,  however,  undergo  clustering,  charge  exchange  and 


ion  atom  Interchange  processes  which  are 

Ng+  +  Ng  +  Ng  -♦  N4+  +  Ng  (38) 

N+  +  Ng  +  Ng  -»  Na+  +  Ng  (39) 

0g+  +  Og  +  Og  ->  04+  +  Og  (40) 

0g+  +  Ng  +  Ng  -»  02+  •  Ng  +  Ng  (41) 

N2+  +  Og  ->  0g+  +  N2  (42) 

N+  +  Og  -*  0g+  +  N  (43) 

0+  +  02  -♦  02+  +  0  (44) 

0+  +  Ng  -4  NO+  +  N  (45) 

N+  +  Og  -*  N0+  +  0  (46) 
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The  clustered  Ions  will  recombine  with  the  free  electrons  via  dissociative 
recombinations  e.g. 


N4*  +  e  -*  N2  +  N2 

(Vo 

N3*  +  a  «♦  N  +  N 

(48) 

02+  ‘  N2  +  e  -*  N2  +  Og 

(49) 

O4  +  e  *♦  Og  +  Og 

(50) 

and  by  three-body  recombinations  (see  Reactions  32-35) 
These  clusters  1  however,  undergo  charge  exchange, 
switching  processes 

rearrangement  and 

N4*  +  0a  ->  Og*  +  2  Ng 

(5D 

/  NO*  +  Ng  +  0 

N3  +  Og  4  V 

'  NO*  +  Ng 

(52) 

Og*  •  N2  +  02  04*  +  Ng 

(53) 

It  Is  obvious  from  reactions  1*2 -46,  51  and  52  that  there  are  various 
charge  exchange  and  Ion  atom  rearrangements  which  proceed  to  produce  Og  and 
No"*",  where  the  last  Ion  recombines  dissoclatively  with  the  plasma  electrons 
(see  Reaction  25)  and  by  various  three-body  recombinations  (see  Raactlon  32- 
35) •  Finally,  the  mutual  neutralization  of  the  positive  and  the  negative  Ions 
can  proceed  via  two  and  three  body  processes.  The  former  Is  represented  by 
the  following  reactions 

0g“  +  A*,  Ag*,  A  *  «*  neutral  products  05*0 

n 

o"  +  A+,  Ag*,  A  *  neutral  products  (55) 

n 

■f*  *f 

where  A  ,  A2  and  Ar  are  atomic,  molecular  and  cluster  Ions,  respectively. 

For  the  three-body  neutralization  add  M  to  both  sides  of  equations  54  and  55. 


14 


However,  when  air  la  moist  and  contains  some  amount  of  HsO,  the 

deionization  becomes  more  complicated  as  discussed  In  the  next  section. 

3.  DEIONIZATION  IN  AIR  WITH  WATER  VAPOR 

The  presence  of  water  vapor  in  air  adds  additional  complication  to  the 

deionization  of  the  disturbed  air*  Various  positive  clustered  ions  are  formed 

with  the  terminal  ion  being  a  hydrated  hydronlum,  H3O  •  (HgO)  .  These  ions 

have  been  observed  by  Narcisl  and  Bailey29  and  the  reactions  leading  to  the 

*4* 

generation  of  these  hydrates  starts  with  Os  as  a  precursor  ion  and  is  well 
understood?7*28  These  hydrates  are  also  generated  starting  with  the  hydration 
of  NO**”, 

The  pertinent  reactions  which  lead  to  the  generation  of  various  hydrated 

4*  *4*  *4*  4* 

ions  can  be  summed  as  follows.  The  positive  ions  N  ,  N<a  ,  N3  ,  N4  and  0 
react  with  HgO  according  to  the  following  reactions. 


N*  +  Ha0  -»  HbO*  +  N 

(56) 

Na+  +  Ha0  -»  H80+  +  Na 

(57) 

Na+  +  Ha0  4  NaH+  +  OH 

(58) 

Na+  +  H^  4  HjgNO+  +  Na 

(59) 

N4+  +  H20  4  HgO"*"  +  2Na 

(60) 

0+  4-  Ha0  4  HbO"*”  +  0 

(61) 

4. 

Reactions  (56)  through  (6l)  obviously  produce  predominantly  HjjO  which  in 
turn  undergoes  the  following  reactions 

(62) 
(63) 
(610 

Reaction  generates  the  hydronlum,  H30~,  which  in  turn  generates  the 
hydrated  hydronlum  according  to 


HaO+  +  02  0a+  +  H20 

HsO+  +  e  -»  OH  +  H 
Ha0+  +  Ha0  *  H*0+  +  OH 
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(65) 


H30'f  +  H2C  +  M  -»  H30+  •  HsO  +  M 


The  following  reactions,  an  association  and  two  switching,  result  in  the 
production  of  02+  •  Hj>0 


02+  +  H2O  +  M  -»  02+  •  U20  +  M 

(66) 

04+  +  H20  4  02+  •  H20  +  M 

(67) 

02+  •  N2  +  H20  -»  0a+*  HaO  +  N2 

(68) 

On  the  other  hand, 
turn  produces  H30+ 

4» 

02  •  HsO  reacts  with  HsO  to  produce  H30  •  OH  which 

•  H20  through  a  switching  reaction  with  H&0.  This 

in 

indicates  two  routes  for  the  production  of  H30+  •  H&0,  one  Initiated  by 

"f*  4*  "t*  4* 

H30  ,  end  the  other  by  02  *  H20  where  0£  and  04  are  precursor  Ions.  Other 

hydrates  are  generated  according  to 

H30+  •  (H20)n  +  HaO  +  M  4  H30+  •  (Ha0)ttfl  +  M  (69) 

All  these  new  clusters  recombine  with  the  free  electrons  through  the 
dissociative  recombination  and  mutual  neutralization  process  to  produce 
neutral  products* 

As  for  the  Impact  of  water  vapor  on  negative  Ions,  the  attachment  of  the 
free  electrons  to  0S  proceeds  as  In  Equation  (36)  with  M  -  HsO. 

The  negative  Ions  also  generate  their  own  ion  clusters  In  the  form  of 
NO3  ,  C03  ,  etc.,  which  are  subsequently  converted  into  hydrated  negative 
ions. 

REACTION  RATES  AND  COEFFICIENTS 

To  calculate  the  air  plssma  conductivity,  a  time  dependent  model  for  the 
electron  energy  deposition  coupled  to  a  detailed  air  chemistry  code  is 
required  .  The  electron  energy  deposition,  however,  requires  a  large  set  of 
Inelastic  cross  sections  for  electron  collision  with  the  air  species.  A  set 
of  cross  sections  appropriate  to  this  problem  is  reported  elsewhere.1* 

A  detailed  air  chemistry  code,  however,  requires  the  appropriate  reaction 
ratoa  which  determine  the  time  history  of  each  species.  Furthermore,  one 
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requires  the  calculation  o£  the  plasma  electron  temperature ,  the  g&a  kinetic 
temperature,  assuming  the  ions  and  neutrals  to  have  the  some  temperature,  and 
vibrational  temperature  o£  the  molecules.  This  is  essential  because  the 
various  reaction  rates  presented  in  the  preceding  sections  depend  on  these 
temperatures. 

In  this  section  we  review  the  pertinent  reaction  rates,  coefficients  and 
their  temperature  dependence  in  order  to  delineate  their  utility  for  the 
modeling  of  the  air  plasma.  It  is  obvious  that  many  of  these  coefficients 
are  known  over  a  limited  range  of  temperatures  and  their  extension  to  higher 
temperatures  may  not  be  accurate. 

4.1  THE  DISSOCIATIVE  RECOMBItVriON  OF  AIR  IONS 

^  iji  .j* 

The  dissociative  recombination  rate  coefficients  of  N2  ,  Og  and  NO 
have  been  measured  under  several  distinct  and  different  experimental  methods. 
These  are,  the  afterglow,30’ 31  the  ion  storage32  and  the  merged  electron-ion 
beam  experiments,33  These  different  techniques  yield  results  in  good 
agreement  for  the  dissociative  recombination  of  02*.  They  predict  30»38»33 
an  electron  temperature  dependant  rate  coefficient  which  varies  as  T^"0,8  for 
an  electron  temperature  from  0,1  to  1,0  eV. 

For  N0+,  however,  the  agreement  on  the  measured  31"33  dissociative 
recombination  rate  coefficient  and  its  dependence  on  the  electron  temperature 
is  not  satisfactory.  For  example,  the  afterglow  experiment31  yields  an 
electron  temperature  dependence  of  x^-0’37  while  the  ion  storage  technique32 
gives  a  temperature  dependence  of  t^"0*33,  The  merged  electron  ion  beam 
experiment33  predicts  a  temperature  dependence  of  T^"0’5,  There  la  a 
reasonable  agreement  between  the  results  of  the  ion  storage  technique32  and 
the  merged  electron- ion  beam  experiment?3  However,  these  results  differ  by  a 
factor  of  2  to  3  from  thone  of  the  afterglow  measurements31  in  the  electron 
temperature  range  of  0. 1-1.0  eV. 

For  N2  ,  the  afterglow30  and  the  merged  electron-ion  beam  method33  yield, 
for  the  dissociative  recombination  rata  coefficient,  a  temperature  dependence 
of  t#"0*39  and  T#"°*s,  respectively.  The  disagreement  is  30#  in  the 
electron  temperature  range  of  0.034-1.0  eV, 

The  disagreement  in  the  dissociative  recombination  coefficients 
summarized  in  Table  1,  may  bo  due  to  the  ions  being  vlbrationally  excited31 
compared  to  those  being  in  the  ground  state.33 


Table  1  The  Dissociative  Recombination  Rate 
Coefficients  of  Ng1",  0a+  and  N0+ 


Na+  +  e  •*  N  +  N 


+  e 


0+0 


+  e 


N  +  0 


(T  in  eV) 

6 

4.3  X  10"8  T 

Q 

5.5  X  10"a  T  “°  * s 

Q 

1.5  X  10"8  T  ■°*7,  T,  £  0.1 

6  t 

2.1  X  10"8  T  "°*5,  T  >  0.1 
e  *  e 

2  X  10“8  T  ■°‘8a 

ft 

9.8  X  10“8  T  "0,a7 
6 

3.6  X  10"8  T  ' “°*5 

ft 


Reference 

30 

33 

30 


32 
31 

33 


4.2  DISSOCIATIVE  RECOMBINATION  RATE  COEFFICIENTS  OF  AIR  ION  CLUSTERS 

Ofa  i.  ^,4.,.  i— .u.) 

The  rate  coefficient  for  the  dissociative  recombination  of  N4  has  been 
measured  at  room  temperature  (T  ■  T  )  by  Kasner  and  Biondi.34  There  Is 
no  temperature  dependent  measurement  for  the  recombination  coefficient , 
however,  the  temperature  dependence  has  been  Inferred  from  the  measurement 
of  an  effective  recombination  coefficient35  for  a  pure  Ng  plasma.  This 
effective  recombination,  measured  in  the  temperature  range  of  295  *  694° K 
with  electrons,  ions  and  the  neutrals  having  the  same  temperature,  yields  a 
T a"1*07  dependence  for  the  recombination.  Since  the  measurement35  was  made 

at  nitrogen  pressures  of  10  -  15  Torr,  the  effective  recombination  coefficient 

■f  4* 

is  in  reality  close  to  that  for  the  recombination  of  N4  ,  because  N4  will  be 

*4*  4* 

the  dominant  ion.  However,  other  ions  such  as  Ng  and  N3  are  also  present 
albeit  at  quantities  much  less  than  that  of  N4+.  It  is  of  Interest  to  note, 
however,  that  the  room  temperature  value  for  the  effective  recombination  rate 
coefficient35  is  in  good  agreement  with  the  dissociative  r <  combination  of 
as  measured  by  Kasner  and  Blondl.34 
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There  is  no  measurement  for  the  dissociative  recombination  of  N3+.  It 

is  expected,  however,  that  the  rate  coefficient  for  its  dissociative 

*4*  + 

recombination  to  fall  somewhere  between  those  of  N2  and  N4  .  An  estimated 
value  and  its  temperature  dependence  of  T  “1,1±0,s  has  been  proposed  by 
Biondi. 30,37 

As  for  04+  there  exist  two  measurements  for  its  dissociative 

recombination,  one  at  a  temperature  of  205°  X  (ref.  36)  and  the  other  at 

180°K  (ref.  39).  The  extrapolation  of  this  data  to  higher  electron 

temperature  is  assumed  to  follow  a  T  “1'°  dependence  in  similarity  to  N4+. 

8 

The  dissociative  recombination  rate  coefficients  for  the  air  ions  are 
given  in  Table  2,  where  the  coefficient  for  02+  *  Na  is  only  an  estimate40 


Table  II  The  Dissociative  Recombination  Rate 
Coefficients  of  Air  Ion  Clusters 


Reaction 

Coefficient 

Reference 

N4+  +  e 

■4  Na  +  N2 

(T  in  eV) 

2  X  10"a,  T#  -  0,025  eV 

34 

3.4  X  10“8 

1 

(-7-  )1*1 

35 

e 

N3+  +  e 

■4  N  +  Na 

1.75  x  10"8 

<T>  11 0,5 

3  6, 

37 

e 

04+  +  e 

•4  oa  0a 

2  X  10“®, 

T  -  0.017  eV 

e 

38 

3.4  X  10“a 

(-)1*0 

3  6, 

37 

e 

1.8  X  10“6, 

T  -  0.015  «V 

8 

39 

02+  •  Na 

+  e  «♦  0a  +  N2 

3.7  x  10“8 

40 

e 

CLUSTERS 

The  dissociative  recombination  rata  coefficient  of  hydronlum,  H30  , 
has  been  measured  by  Leu,  et  al41  at  a  temperature  of  540°  K  (electrons,  ions 
and  neutrals  were  in  thermal  equilibrium).  Heppner,  et  al,42  have  measured 
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Che  cross  section  for  the  dissociative  recombination  in  the  electron  energy 

range  of  0.06  *  1.15  eV  where  the  gas  temperature  was  kept  at  kOO°K. 

Averaging  this  cross  section  over  a  Maxwellian  electron  velocity  distribution, 

Heppner,  et  al42  obtain  a  rate  coefficient  which  follows  a  temperature 

dependence  of  T^-1'®  between  0.025  and  1  eV  and  is  in  good  agreement  with  the 

rate  proposed  by  Biondi.38  However,  the  temperature  dependence  follows  Te-1‘2 

above  T  ■  2eV  and  reaches  a  T  -1*43  dependence  at  T  ■  lOeV. 
e  e  e 

In  addition  to  the  above  measurements,  41,4a  the  dissociative 
recombination  of  H3O*  has  been  measured  extensively  in  flames,  shock  tubes, 
and  afterflow  experiments.  Two  previous  measurements43’44  in  the  temperature 
range  of  (2000  -  2500°  K)  and  a  measurement  45  at  1000° K  are  in  good  agreement 
with  results  of  Heppner,  et  al?s  Other  high  temperature  measurements  of  the 
rate  coefficient,  with  T  -  T.  -  T  yield  reasonable  values,  in  the 

6  1  d 

temperature  range  of  2000  -  2500° K,  compared  to  measurements  of  Heppner, 

et  alt®  However,  the  measurement  of  Ogram,  et  al48  yields  a  dependence  of 

T^-o.s  £or  t^e  rate  coefficient  in  the  temperature  range  of  0.075  to  0.77eV. 

This  temperature  dependence  is  similar  tc  that  obtained  using  the  cross 

section  measurement  (T  s  0.08eV)  of  McGowan  et  al47  However,  they  differ 

e 

by  more  than  a  factor  of  2  in  magnitude.  The  summary  of  these  rate 
coefficient  measurements  is  given  in  Table  III. 

4* 

There  is  no  measurement  for  the  dissociative  recombination  of  H2O  in 

afterglow  type  experiment.  This  may  primarily  be  due  to  the  fast  switching48 
4.  4. 

reaction  of  H2O  with  H2O  generating  H3O  .  However,  merging  beam  technique 
predicts  ~  1/^Te  dependence  for  the  recombination  coefficient  for  electron 
temperatures  of  Te  £  0.1  eV. 

The  dissociative  recombination  of  Hs0+  *  (HgO)  has  been  measured41  at 
4-15  and  54-0°K  while  for  HsO  •  (HjgO)a  the  measurements  are  made  at  300,  415 
and  5^0° K*  Those  temperatures  are  thermal  equilibrium  values,  where  T  -  T  . 

6  A 

From  these  measurements  the  temperature  dependence  for  these  dissociative 

recombinations  are  deduced  and  are  given  in  Table  II.  It  should  be  noted, 

4* 

that  the  dissociative  recombination  of  the  hydrated  hydronium,  H3O  •  (Hz0)n 
have  a  weak  or  no  dependence  on  the  electron  temperature,  especially  for 
n  *  2. 

Few  other  clusters  of  interest  e.g.  O2  *  HaO,  and  H3O  *  OH.  have  no 
measured  or  calculated  values  for  their  recombination  coefficients.  However, 
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estimates  can  be  given,  and  are  presented  in  Table  III. 


Table  III 

The  Dissociative  Recombination  of  H3O*  and  Some  Pertinent  Clusters 

Reaction  Coefficient  (cm3/sec)  Reference 

(Te  in  eV) 

H30+  +  e  -»  H  +  H2O  3.2  x  10"7  T  "1*°,  0.025  s  T  s  0.S6  36 

e  e 

1.0  x  10“f  T  -  0.046  41 

e 

2.56  x  10“®  CjT-)1*2,  0.26  s  Te  <:  2.15  42 

e 

3.05  x  10~e  fc}-)1*43,  2.15  <  T  <  8.6  42 

e 

6.0  x  10-7  Te  -  0.086  45 

2.0  x  10“7,  T  -  .172  43 

e 

1.1  x  10*a  T  ”0#s,  T  £  0.08  47 

1.18  x  10'7  Te*°*s,  Tft  s  0.77  46 

Ha0+  +  e  -4  0  +  Ha  1.2  x  10“8  T  s  0.1  eV  47 

e 

H30+  •  HaO  +  e  -♦  2HaO  +  H  1.6  x  10"®  (~-)0*15  36,  37 

Ae 

H30+  *  (HaO) a  +  e  -»  3HaO  +  H  4.2  x  10“®  (~-)0,oa  36,  37 

e 

02+  •  HaO  +  e  -4  0a  +  HgO  7.2  x  10“7  (;jr)0,s  37 

e 

H30+  •  OH  +  e  4  H20  +  H  +  OH  9.5  x  10"7  (^-)0*2  37 

e 

4.4  POSITIVE  ION-NEUTRAL  ASSOCIATION 

The  three-body  association  rate  coefficient  of  Na+  with  2N2  to  form 
N4*  has  been  measured48"59  at  room  temperature,  300°K  and  as  a  function  of 
E/P  (E  is  the  electric  field  and  P  is  the  pressure) .  The  most  recent 
measurement93  yields  a  rate  coefficient  of  5*5  x  10" 89  cm3/ sec.  The 
temperature  dependence  can  be  inferred  from  the  measurements  of  Good,  et  als 
which  ware  made  at  300° K  and  380° K.  An  80^K  rise  in  the  gas  temperature 
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reduces  the  association  rate  coefficient  by  a  factor  of  2.5  indicating  a 
T  -3.9  dependence.  A  most  recent  measurement61  of  the  association  coefficient 

A 

with  the  third  body  being  He  gives  a  value  of  1.1  x  10"29  cm8/ sec.  If  one 
normalizes  this  value  for  N2  as  a  third  body  by  multiplying  it  by  an  empirical 
factor  of  4  determined  by  Niles61  a  value  of  6  x  10-29  cm6/sec  is  obtained. 

A  suggested37  temperature  dependence  of  T  "1,°  has  been  given  with  an 

A 

uncertainty  in  the  temperature  exponent  of  (+1  to  -2.0). 

The  three-body  association  rate  coefficient  of  N  with  2N2  to  form  N3 
has  been  measured  57*s9»a3  at  room  temperature,  with  the  recent  value59  being 
2.7  x  10"29  cm6/sec  compared  to  5,0  x  10*29  cm6/sec  (ref.  57)  and  1.8  x  10 “29 
cm6/sec  (ref.  63).  A  recent  measurement61  of  the  association  coefficient 
with  the  third  body  being  He  gives  a  room  temperature  value  of  5*2  x  lO**30 
cm6/ sec.  If  on  utilizes  the  empirical  factor62  of  ^  to  convert  the 
coefficient  from  He  to  nitrogen  as  a  third  body,  one  obtains  a  value  of 
2.1  x  10“  29  cm8/ sec,  suggesting  that  the  recent  measurement59  to  be  a  best 
value.  There  exists,  however,  no  measurement  at  temperatures  above  300° K. 
Thus,  the  temperaturedependence  for  the  association  of  N+  to  form  N3+  is  only 
estimated.37 

The  three-body  association  rate  coefficient  of  O2  with  202  to  form  O4 
has  been  measured64  at  room  temperature  and  as  a  function65  of  gas 
temperature  (between  Q0°K  and  300°K).  The  room  temperature  rate  coefficient 
is  2.5  x  10_3°  cm6/ sec  and  the  temperature  dependence  inferred  from  this 

limited  range  follows  T  “3*2.  The  association  of  02+  with  2N2  to  form 
02  •  N«  has  an  estimated3 '  rate  coefficient. 

The  association  of  02+  with  HsO  to  form  02  •  H^  with  the  third  body 

being  N2  and  02  has  been  measured64' 66  and  the  corresnonding  rate  coefficients 
at  T  ■  300°K  are  2.8  x  10”29  cm6/sec  and  2.3  x  10"29  cra6/sec,  respectively. 

a 

The  temperature  dependence  of  this  coefficient  is  estimated.37 

The  three-body  association  of  hydronium  to  form  the  first  hydrated 
hydronium  has  been  measured63*65  with  the  third  body  being  N2  and  02.  The 
room  temperature  coefficients  are  3.^  x  10”27  cm6/sec  and  3.7  x  10”27  cms/sec 
for  M  ■  N2  and  M  -  02,  respectively.  However,  the  temperature  dependence  of 
the  association  is  estimated.37 

The  three-body  association  of  ll30+  •  H^  with  HgO  to  form  H30+  •  (Ha0)2 
has  also  been  measured57  with  the  third-body  being  N2  and  02.  The  room 
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temperature  coefficients  are  2.3  x  10"®7  cma/«ec  and  2.0  x  10-arr  cm*Vsec  for 
M  •  Ng  and  M  •  0a,  respectively.  The  temperature  dependence,  however  Is 
estimated.37  The  rate  coefficients  described  in  this  section  are  summarised 
in  Table  IV. 


Table  IV 

Positive  Ion-Neutral  Association  Coefficients 


Reaction 

Coefficient  (cmg/sec^ 

Miaiq 

Ng+  +  2N2  -4  N4+  +  Ng 

5.5  x  10"C;  t^  «  0.025 

59 

8.0  x  10"®*,  T^  -  0,025 

57 

3.2  x  10“2®  T^  -  0.031 

57 

6.0  x  r"30,  T  -  0.025 

Cl 

6o 

1.7  x  10“ M  »  He 

6i 

i  a. s 

3.1  x  l0“aa,  (-j-) 

57,  59 

i*  +  i'rt 

1.3*10-*,  (-f) 

57 

N+  +  2N2  -♦  N3+  +  N2 

2.7  x  10“®*,  Ta  -  0.025 

59 

1.8  x  10“5*  ,  Ta  -  0.025 

63 

5.2  x  10" ^  M  -  He 

6l 

0g+  +  202  -4  04+  +  Og 

2.1  x  10  "2®  ,  M  -  N2 

2.5  x  10"30,  T^  ■  0.025 

see  text 

64 

i  3.2 

1.94  x  io"aB,  m 

X 

35 

i®  3.2 

3.8  x  10"35,  (“)  ,  M  -  N2 

ti 

37 

02+  +  M  +  HgO  <4  02+  *  H20  +  M 

2.8  x  ■10“4d  ,  M  -  N2  -  0.025 

64,  66 

2.3  x  1O"20,  M  -  02  x  -  0.025 

a 

64 

1.75  x  10-31  M  -  N2 

37 
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Table  IV  continued 


Reaction  Coefficient  Reference 

1.43  x  10“31  (-f-)8*0,  M  »  02  37 

a 

0/  +  N2  +  N2  ■+  0a+  *  N2  +  N2  5.6  x  10"34  (T  )“2,°  37 

a 

H30+  +  HjjQ  +  M  -»  H30+  .  HaO  +  M  3.4x  10"27  T  -  0.025,  M  -  Na  57 

A 

3.7  x  10”87,  T  -  0.025,  M  -  02  57 

ft 

2.1  x  10"30  (T  'T2*0  ,  H  -  Na  37 

& 

1.4  x  10“33  (T  )“4,°,  M  -  02  37 

A 

H30+  •  HaO  +  H20  +  M  h  2.3  x  lO*2^  M  -  N2,  -  0.025  57 

H30+  •  (HaO)  4-  M 

2.0  x  I0mzrt  M  -  Os.  Ta  -  0.025  57 

1.4  x  lO*30  (y-)2  »  M  -  Na  37 

a 

7.8  x  10“33  37 

4.5  CHARGE  EXCHANGE.  SWITCHING  AND  REARRANGEMENT  REACTIONS 

The  reaction  of  N+  with  02  has  three  outgoing  channels,  i.e. 


+  02  -» 

02+  +  N 

(a) 

-♦ 

N0+  +  0 

(b) 

-♦ 

0V  +  NO 

(c) 

The  overall  reaction  rate  has  been  measured  at  room  temperature  by  McFarland, 
et  al37  Smith,  et  al,31  and  McCrumb  and  Warnock?9  The  corresponding 
measurements  yield  5*5  x  10-1°  6.1  x  lO"10  and  7.5  x  10“10  in  units  of  cm3/sec, 
respectively.  The  branching  ratios  have  also  been  measured  with  agreement 
within  10#  for  reactions  a  And  b  and  30#  for  reaction  c.  The  measured 
temperature  dependence  of  the  overall  rate  coefficient  yields3 1  no  apparent 
dependence  on  temperature  up  to  T^  «  0.39eV.  However,  above  T^  -  0.39eV 
the  temperature  dependence  follows  T^°*s7  .  Thus,  on?  can  utilise  the  overall 
reaction  rate  and  the  branching  ratio  as  measured  by  Smith,  et  alsl,  couple! 
to  the  temperature  dependence  obtained  by  McFarland,  et  al67  to  obtain  the 
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the  best  x  ,a  for  these  reactions. 

Tiie  rate  coefficient  for  the  charge  exchange  of  Ns**  with  Os  has  bean 
measured  at  room  temperature81* 87  and  as  a  faction  of  the  ion  temperature.87 
The  room  temperature  values  of  reference  6l  and  6j  are  5,1  x  10"11  cm3/sec 
and  5.0  x  10”11  cm3/  ec,  respectively. 

The  reaction  of  N3+  with  Os  resulting  in  two  channels  which  produce  N0+ 
and  N02+,  have  been  measured  at  room  temperature. 59,81,68  The  temperature 

dependence  of  the  coefficient  for  the  overall  reaction  has  also  been  measured87 
in  the  temperature  range  of  0.025eV  to  0.075eV.  The  overall  room  temperature 
values  are  7*7  x  10"11  cm3/sec,  5,1  x  IQ"11  cm'3/sec  and  6  x  10"11  cm3/sec 
as  obtained  by  references  59 ,  6 1  and  68 ,  respectively. 

The  rate  coefficient  for  the  charge  exchange  of  N4+  with  02  has  been 
measured  at  room  temperature59 *81  and  as  a  function69  of  the  ion  energy  in  the 
range  of  0.0^4-  to  0.2eV  (T^  —  0.025  -  0.132eV).  The  room  temperature  values 

of  the  coefficient  are  k  x  IO"30  cm3/sec,  2.5  x  10“10  and  2.5  x  10"10  cm3/sec, 
ae  obtained  by  references  5*»  6l  and  69,  respectively.  No  temperature  depen¬ 
dence  has  been  observed  for  the  ion  temperature  from  0.025  eV  to  0.0132  eV. 

The  rate  coefficient  for  the  charge  exchange  of  0  with  0a  have  been 

measured  at  room  temperature81’87*70  and  as  a  faction  of  the  ion 
temperature.67*70  The  room  temperature  values  are  1.9  x  10"11  cm3/sec. 

2.0  x  10”11  cm3/ sec  and  2,0  x  10"1-1,  cm3/sec  as  obtained  in  references  6l,  67 
and  70 »  respectively. 

•f* 

The  rate  coefficients  for  the  charge  exchange  of  N  ,  N4  and  0  with 
H20  have  been  measured01  recently,  where  the  corresponding  room  temperature 
values  are  2.8  x  10"9  cm3/sec,  3.0  x  10“9  cm3/ sec  and  3.2  x  10  9  cm3/sec, 
respectively.  The  charge  exchai^e  of  Ns  with  HsO,  however,  has  two  channels0"1- 
resulting  in  H;20+  and  N2H+  with  a  branching  ratio  of  0.82  and  0.18, 
respectively.  As  for  N3+,  its  reaction  with  yields81  HsNO"*". 

The  rate  coefficient  for  the  ’•eaction  of  0+  with  N2  to  form  N0+  has  been 
measured  recently  and  the  room  temperature01  coefficient  is  1.2  x  10" 12 
cm3/ sec.  This  is  in  good  agreement  with  previous  measurement':71  of  the  rate 
coefficient.  However,  the  reaction  0+  with  depends  on  both  the 
vibrational  temperature  of  the  N2  and  on  the  kinetic87*71  temperature  of  the 
ion. 

The  rate  coefficients  for  the  preceding  reactions  are  summarized  in 


Table  V  along  with  few  pertinent  switching  and  rearrangement  processes 


Table  V 

Charge  Exchange ,  Switching  and  Rearrangement  Reaction  Rate  Coefficients 


g&as.tioa 

N*  +  02  «♦  0*+  +  N 


N*  +  02  N0+  +  0 


N*  +  0a  -t  0+  NO 


N2*  +  02  -*  02+  +  N2 


N3*  +  02  5  NO*  +  N2  +  0 


i  no2*  +  n2 


<Ta  in  eV) 

JAx  1Q"“  T 

ft 

3.1  x  10"10  T  -  0.025 


Reference 


3.1  x  10"“  T  -  0.025 

ft 

2.6  x  10"“,  T  -  0.025 


6.7  x  10*11,  T 

< 

3.7  x  10"11,  T 


5.1  x  10"11,  T#  -  0.025 


0.025 

59 

0.025 

6l 

0.025  - 

0.39 

37.67 

57  ,  >  0.39 

37,67 

0.025 

59 

0.025 

61 

0.025  - 

0.39 

37,67 

57,  use 

Is  made 

of  data 

from  refs. 

60  and  66 

0.025 

59 

0.025 

61 

0.025 

67 

0.025 

61 

“0  •  0  »r 

*  T* 

£  0.3 

67 

A 

w 

H 

* 

0.3 

67 

0.025 

59 

0.025 

61 

26 


nirtti  i 


Referenc 


Table  V  continued 

mum  coolant 

Q,8  x  10'12  (T/0*52,  Ta  -  0.025  - 

0.0T5 

a,  0.43,  b  -  .57 
a  >  .  5 .  b  ■  .5 
a  a  .7,  b  -  .3 


N4+  +  02  -♦  02+  +  2N2 

4.0  x  10- 10 

2.5  x  10“10 

2.5  x  10'10  v 

o+  +  oa  -*  °2+  +  0 

1.9  x  10’11 

2.0  x  10"11 

2.0  x  10"11 

4.6  X  10“la  (Ta)“0'4  Ta  s  .155 

67 

1.0  x  10“10  (Ta>1,2  Ta  >  0.155 

67 

N+  +  H20  -♦  H20+  +  N 

2.8  x  10"9,  Tfl  -  0.025 

N4+  +  H«Q  -♦  H20+  +  2N2 

3.0  x  10“9,  -  0.025 

0+  +  H20  -»  H20+  +  0 

3.2  x  10“ 9,  Ta  -  0.025 

N2+  +  HaO  ^  H20+  +  N2 

2.8  x  10“9,  Ta  -  0.025 

£  N2H+  +  OH 

a  -  0.82 

b  -  0.18 

n3+  +  h^o  -»  h2no+  +  n2 

3.3  x  10“lof  Ta  ■  0.025 

0+  +  N2  •*  N0+  +  N 

1.2  x  10“12  ,  Ta  -  0.025 

Table  V  continued 


Reaction  Coefficient  Reference 

7-5  x  10‘13,  0.025  T  <  .1  6r(  ,70 

3.2  x  10*11  (T^)1,36  >  0.1 


X  . 

a 

04+  +  HaO  -»  02+  ’  HaO  4-  02  1.5 

02  •  HaO  +  H20  3  H30  •  OH  4'  02 

*  H30+  4-  OH  4-  0a 
H30+  *  OH  +  H20  -»  H30+  •  HaO  4-  OH 
Ha0+  4-  H20  -4  H30+  4-  OH  1.8 

KaO+  4-  Oa  ■*  02+  4-  HaO  2.0 

02+  ’  Na  4-  02  •*  04+  4-  N2  1.0 

Qs"  *  Na  4-  H2O  *♦  0a+  •  H2O  4<  Na 


10'9.  X  “  0.025 

64 

*  a 

1.2  x  10"9.  x  “  0.025 

*  m 

64 

a 


a  -  1.0  x  10“ 8 

b  -  2.0  x  10“10  ,  x  -  0.025 

a 


1.4  x  10“  9 .  x  -  0.025 

*  a 

64 

x  10“9  t  T  -  0.025 

37 

x  10‘10 ,  T  -  0.025 

37 

x  10"9,  T#  ■  0.025 

37 

4.0  x  10"9 ,  T  -  0.025 

A 

37 

4.6  OTHER  REACTIONS 

In  this  section  we  discuss  few  other  reactions  whose  rate  coefficients 
are  summarized  in  Table  VI. 

The  three-body  recombination  of  positive  ions  with  the  third  body  being 
an  electron  or  a  neutral  (reactions  32,  33  and  34,  35.  respectively)  have 
not  been  studied  extensively.  The  theoretical  efforts  have  concentrated  on 
hydrogen  plasmas  and  are  generally  valid  for  low  electron  temperatures. 
However,  theoretical  expressions?6  accurate  to  within  a  factor  of  2  -  5  for 
these  recombinations  can  be  given  and  are  listed  in  Table  VI. 

The  mutual  neutralization  coefficients  have  been  measured  mainly36  at 
room  temperature  and  some  over  a  limited  temperature  range.36  Their  values 
and  their  estimated  temperature  dependences  are  listed  in  Table  VI.  For 
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higher  cluster  Ions,  e  recent  measurement78  yields  values  which  ere  lower  by 
s  factor  of  2  compared  to  the  previously  accepted  values.30 

The  three-body  attachment  of  the  free  electrons  to  0g  to  form  0s~  with 
the  third  body  being  Oa,  Ng  and  HaO  have  been  measured73  at  room  temperature. 
However,  the  attachment  rate  coefficient  when  M  ■  Oa  have  also  been  measured 
In  the  electron  temperature  range73  of  300  -  600°K.  As  for  the  M  ■  HgO  the 

temperature  range  is  even  more  limited78  (T  ■  300  -  400° K). 

e 


Table  VI 

Recombination,  Neutralization  and  Attachment  Coefficients 


£aastl?ng 

Coefficient  * 

<T  and  T.  in  eV) 

fafaransa 

e  +  M  +  (A+,  Aa+) 

*  i  a»B 

5.9  x  10‘31 

36,37 

i 

e 

•  (A,Aa)  +  M 

e  +  e  +  (A+ ,  Ag+ ) 

_  i  4.5 

4.3  x  10“ 37  (— ) 

36,37 

i 

e 

-♦  e  +  (A,  Aa) 

0+  +  0“  4  0  +  0 

4.2  x  10“a  [Tt]“0,B 

36,37 

N+  +  0“  4  N  +  0 

4.1x10“°  [ T#] "° * 5 

36,37 

0a+  +0  4  Oa  +  0 

1.6  x  10“°  [ t  ] * 8 

36,37 

4*  M 

0a  +  02  4  Oa  +  Oa 

6.6  x  10“° 

36,37 

N2  +  Oa  4  Oa  +  Na 

2.53  x  10"°  [T#]“0#8 

36,37 

X+  +  Y"  4  X  +  Y 

1.6  x  10"°  [T^]-0’8 

36,37 

X+  +  Y"  4  X  +  Y 

9.5  x  10“  9  [T^]-0’8 

72 

A+  +  B-  +  M  4 

2.9  x  10" 29  [t  ]-8*5 

a 

36,37 

A  +  B  +  M 


X+  +  Y“  +  M  -»  X  +  Y  +  M  1.0  x  10"89  [T#]"8*8  36,37 

*  Two  and  throe  body  coefficients  are  in  units  of  cm3/sec  and  cm°/sec, 
respectively. 
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Table  VI  continued 


Reaction 

e  +  O2  +  Oa  >4  O2  +  Og 
e  +  O2  +  Ng  Oa  +  Oa 
e  +  02  +  HaO  -»  Oa"  +  KaO 
e  +  Oa  -*  0"  +  0 


3.6  x  10‘31  <■£■)  Exp  (-  ^^52  ) 

e  e 

1.0  x  10“31 

1.4  x  10“a'4 

1.5  (T#)“1,°  x  10"w  Exp  (-  ^-) 


Minim 

73 

73 


73 


5.  COMMENTS  ON  THE  REACTION  RATES  AND  SOME  RECOMMENDATIONS 

We  have  presented  in  eeotion  4  the  currant  statue  of  the  reaction  rate 
coefficients  pertinent  to  the  modeling  of  the  ionization  and  deionization 
phenomena  in  air.  These  rate  coefficients,  however,  can  be  categorized  as 
follows: 

a)  Most  of  the  reaction  rate  coefficients  are  measured  at  room 
temperature  (0.025  aV)  only. 

b)  A  small  number  of  the  reactions  have  no  measured  or  calculated  rate 
coefficients,  av.?n  at  room  temperature.  Their  coefficients  are  thus 
estimated,  btfiiac  on  similarity  and  educated  guesses. 

c)  Many  of  the  rvacMon  coefficients  are  measured  as  a  function  of 

the  temperature  0  the  reactant.  However,  the  temperature  range  for 
these  measurement!'  Is  narrow  and  generally  Is  limited  to  few  times 
the  room  temperature. 

d)  Few  of  the  reaction  coefficients  are  measured  over  a  wider 
temperature  range. 

To  model  an  air  plasma  of  moderate  to  high  temperature  [  T  ~  (1-2)  eV, 

Tj  *  0.1  -  0.5  eV] ,  one  must  utilize  the  current  reaction  rates.  In  doing 
so,  however,  one  extends  the  known  low  temperature  dependence  of  the  relevant 
coefficients  to  regions  of  higher  temperatures.  Such  an  extension  is  not 
always  warranted  or  accurate.  Furthermore,  in  cases  where  the  energy 
deposited  in  air  is  high  and  results  in  the  enhancement  of  the  vibrational 
modes  of  air  molecules  and  their  lone,  additional  complications  may  arise.  It 
is  possible  that  many  of  the  pertinent  reactions  will  be  affected  by  the 
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vibrational  tamparatura  In  ona  way  or  anothar.  Unfortunately ,  vary  llttla 
Is  known  or  has  baan  measured  undar  thasa  conditions  axcapt  for  tha  following 
raaction: 

0+  +  N2*  4  NO"*"  +  N 


which  is  anhancad  by  a  factor  of  —  40,  compared  to  tha  raaction  at  room 
tamparatura,  whan  tha  vibrational  tamparatura  is  (0.3  -  0.4  aV) . 

Claarly  thara  is  a  naad  in  axtanding  tha  measurements  of  many  of  thasa 
raaction  coefficients  to  moderate  and  higher  temperatures.  However,  tha 
following  reactions  merit  immediate  attentions. 

K.  Ag.SQC,UT10N  .mc,T10N5: 

Tha  dependence  of  tha  following  asaociation  reactions  on  tha  gas  kinetic 
temperature  is  needed  for  Tft  up  to  0.1  aV, 

Na+  +  Ns  +  Na  4  N*+  +  N2 
02+  +  02  +  02  4  O4"*"  +  02 
02+  +  0a  +  Na  4  02+  •  N2  +  02 

02+  +  HaO  +  M  4  0a+  •  HaO  4-  M 

The  reverse  of  these  reactions,  l.e.  tha  colllslonal  breakups  are  also 
needed  ana  in  the  same  temper# ture  range.  It  is  also  of  immediate  interest 
to  measure  the  association  rates,  whan  thasa  molecules  are  vlbratlonally 
excited  (Ty  up  to  ~  0.4  aV) . 

b.  = 

The  rate  coefficients  for  the  following  reactions  and  their  dependence 
on  the  electron  temperature  is  needed  for  T  up  to  1-2  aV. 

N**  +  a  4  N2  +  N2 

04"*"  +  e  4  02  +  02 

Oa+  *  Na  +  a  4  n2  +o2 
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N3+  +  a  *»  Na  +  N 


HaO+  +  a  *4  Ha  +  0 
H30+  +  a  -4  H2  +  0 
N0+  +  a  -4  N  +  0 

The  dissociative  recombination  of  H30  has  a  large  number  of  measurements 
and  art  carrlad  out  ovar  a  wlda  ranga  of  tha  alactron  and  gaa  tamparaturea 
(aaa  Tabla  III  aaction  L.3).  Tha  measurements  of  Heppner,  at  al,42  at  thla 
Junctura,  praaant  tha  moat  raaaonabla  coafficiant.  Howavar,  a  racant 
measurement48  of  tha  coafficiant  diffara  from  Heppner,  at  al,4a  both  in 
magnituda  and  tha  form  of  tha  tamparatura  dapandanca.  Thua  it  ia  of  lntaraat 
to  dalinaata  and  claar  up  tha  controvaray  among  thaaa  diffarant 
maaauramanta . 42' 48  * 47 

4. 

Tha  diaaociativa  racombination  of  NO  (aaa  Tabla  I  aaction  4.1)  haa  baan 
maaaurad  by  three  diffarant  axparimantal  techniques.31"33  Thay  pradict  tha 
following  tamparatura  dapandanciaa,  t#"0*37,  T#"°‘8a  and  which  ara 

obvloualy  diffarant.  In  addition  to  tha  axparimantal  afforta,  tha 
diaaociativa  racombination  of  NO  haa  alao  baan  inveatigated 
thaoratically. 75-78  Thara  la  aoma  agraamant  on  tha  thaoratical  bahavlor  of 
tha  coafficiant  in  tha  low  and  high  alactron  tamparatura  ragimeo.  Howavar, 
thara  la  no  agraamant  on  tha  magnituda  of  tha  coafficiant  aapacially  on  tha 
dlaaoclation  racombination  of  tha  vibrationally  axcltad  N0+.  Tharafora, 
it  la  of  lntaraat  to  dalinaata  tha  problama  aaaoclatad  with  tha  diaaociativa 
racombination  of  N0+. 

Undar  thla  category  wa  include  the  naad  for  extending  tha  maaauramanta 
of  tha  diaaociativa  racombination  coafflcianta  of  tha  poaitlva  ion  clusters 
to  higher  alactron  tamparaturea  (T#  >  0.0 6  eV) . 

C.  ELECTRIC  FIELD  DETACHMENT  OF  THE  NEGATIVE  IONS 

It  ia  of  lntaraat  to  know  tha  detachment  rataa  of  negative  Iona 
(0",  0a,  clusters)  undar  tha  influence  of  high  electric  fields. 
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